IMS Collections 

Probability and Statistics: Essays in Honor of David A. Freedman 

Vol. 2 (2008) 42-75 

© Institute of Mathematical Statistics, 2008 
DOI: 10.1214/193940307000000383 

Brownian motion on disconnected sets, 
basic hypergeometric functions, and some 
continued fractions of Ramanujan 

Shankar Bhamidi 1 , Steven N. Evans* 1 Ron Peled^ 1 and Peter Ralph* 1 

University of California, Berkeley 

Abstract: Motivated by Levy's characterization of Brownian motion on the 
line, we propose an analogue of Brownian motion that has as its state space 
an arbitrary closed subset of the line that is unbounded above and below: such 
a process will be a martingale, will have the identity function as its quadratic 
variation process, and will be "continuous" in the sense that its sample paths 
don't skip over points. We show that there is a unique such process, which 
turns out to be automatically a reversible Feller-Dynkin Markov process. We 
find its generator, which is a natural generalization of the operator / i— > \ f" ■ 
We then consider the special case where the state space is the self-similar 
set {iq k : k G Z} U {0} for some q > 1. Using the scaling properties of the 
process, we represent the Laplace transforms of various hitting times as certain 
continued fractions that appear in Ramanujan's "lost" notebook and evaluate 
these continued fractions in terms of basic hypergeometric functions (that is, 
q-analogues of classical hypergeometric functions). The process has as a 
regular instantaneous point, and hence its sample paths can be decomposed 
into a Poisson process of excursions from using the associated continuous 
local time. Using the reversibility of the process with respect to the natural 
measure on the state space, we find the entrance laws of the corresponding 
Ito excursion measure and the Laplace exponent of the inverse local time - 
both again in terms of basic hypergeometric functions. By combining these 
ingredients, we obtain explicit formulae for the resolvent of the process. We 
also compute the moments of the process in closed form. Some of our results 
involve g-analogues of classical distributions such as the Poisson distribution 
that have appeared elsewhere in the literature. 



1. Introduction 

Let T be an arbitrary closed subset of K. There is a well-developed theory of dif- 
ferentiation, integration, and differential equations on T (sometimes refered to as 
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the time scale calculus) that simultaneously generalizes the familiar Newtonian cal- 
culus when T = R and the theory of difference operators and difference equations 
when T = Z (as well as the somewhat less familiar theory of (/-differences and q- 
difference equations when T is {q k :fceZ} for some q > 1). The time scale calculus 
is described in [G], where there is also discussion of the application of time scale 
dynamic equations to systems that evolve via a mixture of discrete and continuous 
mechanisms. 

Our first aim in this paper is to investigate a possible analogue of Brownian 
motion with state space an arbitrary closed subset of R. A celebrated theorem 
of Levy says that Brownian motion on R is the unique R- valued stochastic process 
(6)teK+ sucn that: 

(I) £ has continuous sample paths, 
(II) £ is a martingale, 
(III) — i)teR + is a martingale. 

A similar set of properties characterizes continuous time symmetric simple random 
walk on Z with unit jump rate: we just need to replace condition (I) by the analogous 
hypothesis that £ does not skip over points, that is, that all jumps are of size ±1. 
Note that for both R and Z the Markovianity of £ is not assumed and comes as a 
consequence of the hypotheses. 

We show in Section 2 that on an arbitrary T that is unbounded above and below 
there exists a unique (in distribution) cadlag process £ that satisfies conditions (II) 
and (III) plus the appropriate analogue of (I) or the "skip-free" property of simple 
random walk. Namely: 

(F) for states x < y < z in T and times < r < t < oo, if cither £ r = x and £t = z 
or £ r = z and £t = x, then £ s = y for some time s between r and t. 

Moreover, we demonstrate that this process is a reversible Feller-Dynkin Markov 
process with a generator that we explicitly compute. The proof of existence is via 
an explicit construction as a time change of standard Brownian motion. The proof 
of uniqueness (which was suggested to us by Pat Fitzsimmons) relies on a result of 
Chacon and Jamison, as extended by Walsh, that says, informally, if a stochastic 
process has the hitting distributions of a strong Markov process, then it is a time 
change of that Markov process. 

As well as establishing the existence and uniqueness of the Brownian motion on 
T in Section 2, we give its generator, which is a natural analogue of the standard 
Brownian generator / i— > \f" ■ Note that a simple consequence of (II) and (III) 
is that £ has the same covariance structure as Brownian motion on R, that is 
- E B K.]E a &] = s At for all x e T. 

The assumption that the state space T is unbounded above and below is neces- 
sary. To see this, first note that T cannot be bounded above and below, because 
this would imply that if £o = x, then lim^oo E[£ t 2 — t] — — oo ^ x 2 , contradicting 
property (III). Assume now that T is unbounded above and bounded below with 
inf T = a > — oo. Suppose £o = %■ Choose b G T with x < b. Put T = inf{t > : 
£t ^ [a, b)}. Note by the right-continuity of £ and property (I') that £t = b on the 
event {T < oo}. By properties (F) and (III), E[t AT] = E[£ t 2 AT ] - x 2 < a 2 V b 2 - x 2 , 
and so T is indeed almost surely finite. By properties (F) and (II), (£*AT)teK + is a 
bounded martingale with £ tA T = b for t > T almost surely, but this leads to the 
contradiction b — lim^oo E[£ tA <r] = x. The proof that T cannot be bounded above 
and unbounded below is similar. 



11 



5. Bhamidi, S. N. Evans, R. Peled and P. Ralph 



The process £ is constructed as a time-change of standard Brownian motion, 
a class of processes described in Ito and McKean [34], and that has been studied 
variously as "gap diffusions" [44], "quasidiffusions" [7, 46, 47, 48, 49, 50, 55], and 
(one-dimensional) "generalized diffusions" [65, 66, 67, 68]. The process £ that we 
study is a quasidiffusion, so results on quasidiffusions apply in this context - but 
it is a distinguished quasidiffusion among the many possible quasidiffusions taking 
values in T. Quasidiffusions can exhibit behavior considerably different from that 
of £ - for instance, Feller and McKean [22] described a quasidiffusion that has 
all of K as its state space, but spends all its time in Q . These processes (with 
killing and appropriate boundary conditions) were shown by Ldbus [54] , extending 
work by Feller [15, 18, 20] to be the only Markov processes taking values in K whose 
generators are in some sense local, and satisfy a certain maximum principle. Various 
authors [47, 48, 51] have given beautiful spectral representations of quasidiffusions 
using Krcm's theory of strings [13, 19, 41]. 

It is natural to ask about further properties of the Brownian motion on T. In 
the present paper we pursue this matter in a particularly nice special case, when 
T = T q := {±q k : k G Z} U {0} for some q > 1. In this case, the process £ started 
at x has the same distribution as the process (^■£ q 2fc t ) tg K + when £ is started at 
q k x for k G Z. This Brownian-likc scaling property enables us to compute explicitly 
the Laplace transforms of hitting times and the resolvent of £ in terms of certain 
continued fractions that appear in the "lost" notebook of Ramanujan. We can, in 
turn, evaluate these continued fractions in terms of basic hypergeometric functions 
(where, for the sake of the uninitiated reader, we stress that "basic" means that 
such functions are the analogues of the classical hypergeometric functions to some 
"base" - that is, the g-series analogue of those functions). We recall that, in general, 
a g-analogue of a mathematical construct is a family of constructs parameterized 
by q such that each generalizes the known construct and reduces in some sense to 
the known construct in the limit "q — > 1" . This notion ranges from the very simple, 
such as (q n — l)/(q — 1) being the g-analoguc of the positive integer n, through to 
the very deep, such as certain quantum groups (which are not actually groups in the 
usual sense) being the "g-deformations" of appropriate classical groups [9, 35, 40]. 

For a very readable introduction to g-calculus see [36], and for its relation with 
g-series, see the tutorial [45], or the more extensive books [30] or [2]. What we need 
for our purposes is given in Section 11. 

The interplay between g-calculus (that is, g-difference operators, ^-integration, 
and g-difference equations), g-series (particularly basic hypergeometric functions), 
and probability has been explored in a number of settings both theoretical and 
applied. The recent paper [4] studies the connection between g-calculus and the 
exponential functional of a Poisson process 



where N t is the simple homogeneous Poisson counting process on the real line. A 
purely analytic treatment of the distribution of I q using g-calculus is given in [3] . It 
is interesting to note that the same functional seems to have arisen in a number of 
applied probability settings as well, for example, in genetics [10] and in transmis- 
sion control protocols on communication networks [11]. In [42] the Euler and Heine 
distributions, g-analogucs of the Poisson distribution, are studied: distributional 
properties are derived and some statistical applications (such as fitting these dis- 
tributions to data) are explored. These analogues have arisen in contexts as varied 




*dt, q<l, 
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as prior distributions for stopping time strategies when drilling for oil and studies 
of parasite distributions, see the references in [42]. The g-analogue of the Pascal 
distribution has also been studied in the applied context, see [43]. The properties 
of g-analogues of various classical discrete distributions are also surveyed in [52]. 
Both I q and the Euler distribution appear in Section 6, where they come together 
to form the distribution of a hitting time. 

Probabilistic methods have also been used to derive various results from q- 
calculus. A number of identities (including the g-binomial theorem and two of 
Euler's fundamental partition identities) are derived in [59] by considering pro- 
cesses involving Bernoulli trials with variable success probabilities. Several other 
identities (for example, product expansions of g-hypergeometric functions and the 
Rogers-Ramanujan identities) are obtained in [58] using extensions of Blomqvist's 
absorption process. Some properties of g-random mappings are explored in [57]: 
in particular, the limiting probability that a g-random mapping does not have a 
fixed point is expressed via a g-analoguc of the exponential function. Connections 
between g-series and random matrices over a finite field (resp. over a local field 
other than R or C) are investigated in [27, 28, 29] (resp. [1, 14]). 

2. Brownian motion on a general unbounded closed subset of R 
2. 1 . Existence 

Let T be a closed subset of R that is unbounded above and below (that is, bilaterally 
unbounded). We now show existence of a Feller-Dynkin Markov process satisfying 
conditions (F), (II) and (III) by explicitly constructing such a process as a time- 
change of Brownian motion. Let (Bt)teR + be standard Brownian motion on R and 
let be its local time at the point a £ R up to time t > 0. We choose a jointly 
continuous version of £ and we adopt the normalization of local time that makes 
I a family of occupation densities for the Brownian motion; that is, J Q f{B s ) ds — 
J R /(a)£" da for all bounded Borel functions /. Equivalently, for each a the process 
(£")t S R + is the unique continuous non-decreasing process such that (\B t — a\ — 
(■t )tea+ i s a martingale. 

We introduce the following notation from [6]. For a point x £ T set 



If p(x) x say that x is left- scattered, otherwise x is left-dense, and similarly 
if <r{x) ^ x say that x is right scattered, otherwise x is right-dense. Denote by 
T aa ,T s d,T(js and Tad the left and right scattered, left-scattered right-dense, left- 
dense right-scattered and left and right dense subsets of T, respectively 



interval of points in R that arc closer to x than to any other point of T. Thus p is 
the push-forward of m by the m-a.e. well-defined map that takes a point in R to 
the nearest point of T. Note that the support of p is all of T. Define the continuous 
additive functional 



p{x) := sup{y £ T : y < x}, 



<j{x) := inf{y £ T : y > x}. 





and let 9^ be its right continuous inverse, that is, 



0? := M{u : > t}. 
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By the time change of B t with respect to the measure [i ([34], §5, or [GO], III. 21) 
we mean the process £t := B e ». It is easily seen that £ has T as its state space, 
and if Bq = x G T then £o = x also. Moreover, it is not hard to show that £ is a 
Feller-Dynkin Markov process on T. 

We will need the generator of £. For that purpose, we introduce the following 
notation. Write 

Co(T) := {/ : T — ► M : f is continuous on T and tends to at infinity}. 

Define a linear operator Q on Cq(T) as follows. For ieT, set y xr := p(x — r) and 
z x r := o~(x + r). Put 



(G/)(aO 



lim ( Zx ' T _ X f(y x ,r) + X f(Zx,r) - f(x) ) / ((x - y x ,r)(Zx,r - x)) 

riO y Z x r Ux,r %x,r Ux,r / / 

fiVx.r) f(x) f(z x ,r) 



f{yx. x ) - f(x) t /(^x.r) - /(») 

lim 



r Io \ 2/x,r) \%x,r Vx,r) (,%x,r XjyZ x ,r Vx^r) 

on the domain Dom(£/) consisting of those functions / G Co(T) for which the limits 
exist for all x G T and define a function in Cq(T). 

Note that ^ is a natural analogue of the standard Brownian generator / i— > i/" 
and coincides with this latter operator when T = M. Note also that if / is the 
restriction to T of a function that is in Cq(K), then / G Dom((?) and 



(Gf)(x) 



(x-p(x))i>(x)-p(x)) (2;-p(x))(<T(x)-x) + (<r(x)-x)( [ r(x)-p(x)) ' ^ fc 11 ' 
f(p(x))-f(x) f{x) -rpT, 
(x-p(x)Y 2 h x-p(x) ' Xfc 



/'(*) , /WHM 



V(x)-X ' (ct(x)'-x)2 ' ' 



(/,/■ 



Proposition 2.1. TTie £irae change £ of standard Brownian motion B with respect 
to the measure /i is a Feller-Dynkin Markov process on T that satisfies conditions 
(I'), (II) and (III). The generator oft; is (Q, Dom(C?)). 

Proof. We have already noted that £ is a Feller-Dynkin Markov process. Given 
x G T, write F x for the distribution of £ for the initial condition £ = x, and denote 
the corresponding expectation by K x . Under any the property (I') is clear from 
the fact that the support of /i is all of T. Before establishing properties (II) and 
(III) under any P x , we first show that the generator of £ is (Q, Dom(§)). 

Write (H, Dom(7i)) for the generator of £. We begin by showing that (H, 
Dom(W)) = (</,Dom(</)). For .t G T and r > 0, set T xr := vaf{t : d{£ t ,x) > r}. 
By Dynkin's characteristic operator theorem [GO], III, 12.2, / G Dom(H), if and 
only if 

(2.1) „ m ^[/(fc.J]-/M 

1 ' rlO E*[T x>r ] 



exists at every ieT and defines a function in Co(T), in which case this function 
is Hf. 
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Set y Xi r := p(x — r) and z x _ r := a{x + r). Because the support of p is all of T, 
9%, = m£{t e R+ : B t G {y x ,r, z x ,r}} U x , r . 



Thus 



and 



c ^ ~\ Z x .r X 



Consequently, 



Hence it is enough to show for all a; <G T and r > that 
(2.2) E K [T X)r ] = {x-y x>r ){z x , r -x). 

Now 

?V = / ^, p A*(do)= / l a jj xri i{da), P^-a.s., 

and in particular, 

E*[T x , r ] = / E x [£^ r }p(da) 

2(q - gfa.rX^r - g) 
(y x>I .,ic] ^a;,r Dx,r 



fi(da) 



f 2(x — y x r){z x r — a) , , . 
(2.3) + / /'_ ' " f (**) 

(z x ,r — x) [ (a- y x ,r) v{da) 



z x,r Vx,r \ J (y Xir ,x] 



+ (x- y x , r ) / (z X) r - a) p(da) 

J(x,z x , r ) 

But, as we now show, for any points it, v € T, u < v, 

<j{u) — u v — p(v) 



(2.4) / p(da) = v - u - 

f v 2 u 2 a{u) - u v - p(v) 

(2.5) j ^ a p(da) = — - — - u v 

(note the similarity to Lcbcsguc integration up to boundary effects). Substituting 
this into (2.3) gives (2.2) after some algebra. 

Let us prove the identities (2.4) and (2.5). For simplicity, we prove them in the 
special case when T^s H (u, v) = 0. The proof of the general case is similar. Fix 
u, v € T, u < v. Since T is closed, we can write (u, p(v)) \ T as a countable union of 
disjoint (non-empty) open intervals {(a n ,b n ) : n € N}. We note that for any such 
interval (a n ,b n ] C (u, v) and 

f i j \ °(b n )-a n a(a n )-a n a(b n )-b n 
/ M« a ) = o = bn ~ an o + o ■ 
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Summing up over all these intervals the boundary effects cancel telescopically (since 
Td s H (u, v) = 0) and, since a(a n ) = b n , we get 



„ J(a„,b n ] 

dm 



cr(u) — u v — p(v) 
2 2 

cr{u) — u v — p(v) 

U„K,6„] U "' 2 2~ ' 

where again m is Lebesgue measure. Identity (2.4) now follows since (u, v) 
\J n {a n , b n ] U ((it, v) n Tdd) and, by the definition of p, 



p(da) = / p(da) + / p(da) 

(u,v) J(u,v)nr M J\Jja n ,b n ] 

<y(u) — u v — p(v) 
dm . 



<(u,v) ^ - 

To prove identity (2.5), we note similarly that for any n E N 



i j \ l v(bn)—an K a 2 n a(a„)-a„ a(b„)-b r , 



(o n ,6„] 2 2 



so that again 



" ^(a„,6„] „ V 2 

= / 



er(it) — u v — p(v) 
u v 



(t(u) — u v — p(v) 
a dm — u v . 



and (2.5) follows since 



a p(da) = / a p(da) + / a p(da) 

(u,v) J(u,v)m dd J\J n (an,bn] 

a(u) — u v — p(v) 
a dm — u v . 

u,v) 2 2 

By the Markov property of £, in order to show (II) and (III) it suffices to show 
that E x [&] = x and E x [£ t 2 ] = x 2 + 1 for all x e T and t e 1 + . By Dynkin's formula 
[GO], III. 10, for any / € Dom(£) 

( 2 -6) M t :=/(&)- f {Gf)(t.)da 

is a martingale (for each starting point). Note that if we formally apply the ex- 
pression for Qf to f(x) = x (resp. f(x) = x 2 ), then we get Qf{x) = (resp. 
Gf(x) = 1), and this would give properties (II) and (III) if x i— ► x and x i— ► x 2 
belonged to the domain of Q. Unfortunately, this is not the case, so we must resort 
to an approximation argument. 

Fix igT. Given any r > 0, for R > r sufficiently large we have [p(x — r), a(x + 
r)] C (p(x — R), <t(x + R)). For any such pair r, R, there are functions g,h £ Dom(Q) 
such that g(w) = w and h(w) = w 2 for w £ [p(x — R),a(x + R)], and hence 
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Gg{w) = and Qh(w) = 1 for w G [p(x — r), cr(x + r)]. It follows that (£tAr Cir .)teR+ 
and (£ 2 a t x — t A X , a; ,r)teK + are both martingales under P 21 . 
Hence, if < r' < r" , then 

® x kut^„ -e*AT x , r ,) 2 ] =m^AT^ r ,} -m^AT^j 

= E X [t A T x ,r"\ —^f[tA T x>r >] . 

Thus ^tAT x r converges to £ t in L 2 (P a: ) as r — ■> oo, and so 

E*[&]= lim E^ tAT _] = 3 ; 

r — >oo 

and 

as required. □ 

Remark. It is more standard, but we believe less natural (and equivalent) to 
instead regard the generator an operator on continuous functions of R that are 
linear outside the support of /i. Such a generator has the natural interpretation 
2 3u lb ■ These operators appear in relation to diffusion processes in Ito and McKean 
[34], §§5.1-5.3, and were studied further by Feller [15, 16, 17, 18, 19, 20, 21] (where 
the support of a is connected), Lobus [54, 55], and Freiberg [24, 25, 26] (where p, 
is atomless). 



2.2. Uniqueness 

We next establish a uniqueness result that complements the existence result of 
Proposition 2.1. 

We will apply the following result, which is a slight variant of Corollary 3.5 of 
[64] extending results of [8]. We make the assumption that the Markov process 
X is a right process and that the process Y is defined on a space satisfying the 
usual conditions to avoid listing Walsh's assumptions. We also state the result in 
terms of bounded rather than finite stopping times, but this is readily seen to be 
sufficient. The result says, roughly speaking, that if a process has the same state- 
dependent hitting distributions as some strong Markov process, then the process is 
a time-change of that Markov process. (For example, a consequence of the result 
is the celebrated result of Dubins and Schwarz that any continuous martingale is a 
time change of Brownian motion, from which Levy's characterization of Brownian 
motion that we mentioned in the Introduction is an immediate corollary.) 

Theorem 2.1. Let X = (n,J 7 ,!Ft,Xt,8t,'P x ) be a Borel right process with Lusin 
state space E. Assume that the paths of X are cadlag and that X has no traps 
or holding points. Let Y be a cadlag process with state space E that is defined on 
a complete probability space (E,C/,Q) equipped with a filtration (Gt)tew. + satisfying 
the usual conditions. Assume Y = xq for some Xq € E and that almost surely 
the sample paths of Y are not constant over any time interval. Given a Borel set 
B C E , put Sb '■= inf{t > : X t G B} and define the corresponding hitting kernel 
by ttb(x, A) := P x {Xs B € A} for x G E and A C E Borel. Given a bounded 
{Qt)tem + -stopping time T, put r = inf{i >T:Y t (z B}. Suppose that 



Q{Y T G A | Qt} = ttb(Yt, A) 
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for all bounded (Gt)teR+ -stopping times T and all Borel sets A and B. Then there ex- 
ists a perfect continuous additive functional for X with continuous inverse (T't) tg R + 
such that (lr t )teK+ h as the same distribution as (X t )t^R + under P§. 
Proposition 2.2. Let C be a cadldg T -valued process such that Co = z £ T. Suppose 
that C satisfies the counterparts of properties (F), (II) and (III) with £ replaced by 
C- Then C possesses the same distribution as the particular Feller- Dynkin process £ 
of Proposition 2. 1 has under ¥ z . 

Proof. We wish to apply Theorem 2.1. Unfortunately, the process £ has holding 
points unless T = M. We adapt an artifice presented in Remark 1 after Theorem 
3.4 in [64] to circumvent this difficulty. 

Without loss of generality, we may suppose that C is defined on a complete 
probability space (S, G, Q), that this probability space is equipped with a filtration 
(<7t)teR + satisfying the usual conditions, and that (Ct)teR+ an d {Ct ~ 0*eR+ are both 
martingales with respect to (Gt)tem. + - We will use Q[-] to denote expectation with 
respect to the probability measure Q. 

We first show that the sample paths of C do not get trapped forever in any 
state. Given a,b £ T with a < x < b, put R = ml{t > : ( t ^ By the 

counterpart of property (F), CtAR <= b] and hence, by the counterpart of property 
(III), Q[t A R] = Q[C t 2 Ai?] -x 2 <a 2 \J b 2 - x 2 . Thus Q[R] < oo and, in particular, 
Q{i? < oo} = 1. Since this is true for all a and b, it follows that almost surely there 
does not exist a time s £ R+ and a state y £ T such that Q = y for all t > s. 

Let S be a finite (Gt)t£& + stopping time, and put T := mi{t > S : Ct ^ Cs}- It 
follows from the above that T < oo almost surely. Moreover, by the counterparts of 
properties (F) and (II) for ( and the right-continuity of paths, (t £ {p((s), (T (Cs)} 
almost surely with 



Q{Ct = p((s) I Gs} = 

and 

Q{Ct = o-(( s ) | Gs} 



v((s) - p(Cs) 
Cs - p(Cs) 



v«s) - p(Cs) 

on the event {Cs £ ^\^dd}- Thus Ct = Cs almost surely on the event {Cs £ T\T SS } 
and hence, by the counterpart of property (III), S = T almost surely on the event 
{Cs S T\T SS }. On the other hand, it is certainly the case that S < T almost surely 
on the event {Cs G TT SS }. 

We next claim that, conditional on Gs^ the random variable T — S is exponen- 
tially distributed with expectation (Cs — p(Cs))( cr (Cs) ~ Cs) (where the exponential 
distribution with expectation is of course just the point mass at 0). This must 
be so, of course, if C has the same distribution as £, and it is the key to adapting 
Theorem 2.1 to our setting in which the processes involved do have holding points. 
To see the claim, define a function $ : T x T ^ I by 



*(x,y) 



x£T s 



' (y-p(x))(cr(x)-y) 
(x-p(x))(a(x)-x) ' 

0, ieT\T s 



Note that for each fixed x the function ^(x, •) is quadratic. It follows from coun- 
terparts of properties (II) and (III) for C that the process 



M t := l{Cs G T ss } 



tf\{T-S) 

^(Cs,C(S+*)at) + 77 77 xw f , x ~ r 

(Cs - p(Cs))(ct(Cs) - Cs) 



t £ 
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is a martingale with respect to the filtration (£7s , +t)teR+ ■ Note that 

rt 1{T-S>u] 



M t = i{Cs e T ss } 

Hence 



1{T-S>t}+ ( 
Jo 



o (Cs-p(Cs))MCs)-Cs) 



du 



HO e T„( (Q{T - S > , | fa) - 1) = - f l%iI»ag-g>«lM „„, 



o Ks - ("(&))(»(&) - &) 



and so 



l{Cs G T SS }Q{T - S > 1 1 £ s } = l{Cs £ T ss } cxp ( - 



(Cs-p(CsW(cs)-cs)j ' 

as claimed. 

We now apply the device from [64] mentioned above to "embellish" the process 
£ in order to produce a Feller-Dynkin process without traps or holding points. Set 
f := (T ss x M) U ((T \ T ss ) x {0}) C T x M. Put U := mf{t > : £ t ^ £ } and 
C< := t — sup{s < t : £ s 7^ £ t }, with the convention sup0 = 0. That is, C* is the 
"age" of £ in the current state at time t. There is a Feller-Dynkin process (£t, P^'")) 
with state-space T such that under P^-") the process £ has the same distribution 
as the process 

(&,« + *), 0<t<C7, 

(6,Ct), t>tf, 

under P x . 

Fix a Borcl_set BCT, write f B := inf{f G R+ : £ t G 5} for the first hitting 
time of -B by £, and denote by ttb the corresponding hitting kernel. That is, 

7t B ((x,u),A) :=f^{& g eA} 

for (x, a)ef and A a Borel subset of T. It is not hard to see that Tb is finite P^' 11 )- 
almost surely for all (x,u) G T, and hence 7Tb((x, m), ■) is a probability measure 
concentrated on the closure of B for all (x, it) G T. 

Let (-Dt)t<=R + be the analogue of (C t ) tS R + for £. That is, -D t := t— sup{,s < i : Cs 7^ 
Ct}. Given a finite (£/t)teK+ stopping time S, put f := inf{t > S : (( t ,D t ) G B}. 
From what we have shown above, it follows by a straightforward but slightly tedious 
argument that if B is a finite set , then 

(2-7) Q{(Cf,Df) eA\gs} = 7r B ((Cs,D § ),A) 

(in particular, T is finite Q-almost surely). If B is arbitrary, then taking a countable 
dense subset of B and writing it as an increasing union of finite sets shows that 
(2.7) still holds. 

Theorem 2.1 gives that there is a continuous increasing process (Tt)teM + such 
that each T t is a (Gt)t&m. + stopping time, To = 0, and ((Ct, , -Dr t ))teR+ has the same 
distribution as £ under p( 2 -°), (recall that Co = z). In particular, (Cr t )teR+ nas the 
same distribution as £ under P x . Since property (III) holds for £ and its counterpart 
holds for £, we have that (£f. — t)t 6 R + is a martingale and (Cr t — T t ) tg R + is a 
local martingale. Thus (T t — f )teR+ is a continuous local martingale with bounded 
variation, and hence T t = t for all t G R+, as required. □ 
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We note that, by a proof similar to that of Proposition 2.2, one can show any 
cadlag T-valued process with properties (I') and (II) is a time-change of a process 
with the distribution of the process £ constructed in Proposition 2.1. It may be 
necessary to introduce extra randomness in the time-change to convert the holding 
times of the process at points in T ss into exponential random variables, and it may 
also be necessary to introduce extra randomness to "complete" the sample paths 
of the copy of £ - as the original process may "run out of steam" and not require 
an entire sample path of a copy of £ to produce it (the most extreme example is a 
process that stays constant at its starting point). This observation is the analogue 
of the result of Dubins and Schwarz that any continuous martingale on the line is 
a time-change of some Brownian motion. 



2. 3. Reversibility 



Extensions of the following result will hold more generally: under suitable hypothe- 
ses, a time-change of a Markov process that is reversible under some measure will 
be reversible under an appropriate new measure. Since we don't know of a suitable 
general reference, we provide the straightforward proof in our setting where the 
Markov process is Brownian motion. 

Lemma 2.1. The process £ of Proposition 2.1 is reversible with respect to the 
measure fi. In particular, /i is a stationary measure for £. 

Proof. We have to show for all A > and all non-negative Borel functions / and g 
that 



f{x)E x 



5(6) dt 



fj.(dx) 



g(x)E* 



-XI 



/(&) dt 



fj,(dx). 



Now recalling is the inverse of 6* M , 



f{x)W 



e- Xt g(Z t )dt 
f(x)E 
f{x) E 



fx(dx) 



f{x)W 



e- M g(B g »)dt 



fi(dx) 



g(B s )dA» 



[i(dx) 



T JT 



dil 



g(y) v(dy) n(dx) 



It follows from the reversibility of B with respect to Lebesgue measure that for 
any 7 > and any non-negative bounded continuous functions F, G and H , 



F(x)E 
F(x) E x 
G(y)W 
G(y) E y 



£ -hs+\ J T e a s H(a) m (da)) £ 

e - hs+ xJ o H { B u )du) G{Bs)ds 
- e - hs+ xf°H(B u)d u) F{Bs)d - 
g -( 7 s+A f T r s H(a)?n(da)) ^ x 



G(y) m(dy) m(dx) 
m(dx) 
m(dy) 

F(x) m(dx) m(dy). 



Thus (noting that each side is jointly continuous in x and y), 



-(7S+A j^i a s H(a)m(da)) 



dil 



-(7s+A f T £«H(a)m(da)) ,„ x 
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for all i,j£R. 

Writing [i as the vague limit of a sequence of Radon measures that have bounded 
density with respect m and applying dominated convergence gives 



-(~/ S +\ J T e a sP ,(da)) 



Ml 



D -(7S+A f^e" fj,(da)) ,„ x 



for all x,y €M.. Hence, by monotone convergence, 



dil 



for all x, y S R. This suffices to establish the result. 



□ 



3. Hitting times of bilateral birth-and-death processes 



In order to compute certain hitting time distributions for £, we recall and develop 
some of the connections between Laplace transforms of hitting times for a birth-and- 
death process and continued fractions. See Section 12 for some relevant background 
and notation for continued fractions. 

The connection between birth-and-death processes and continued fractions has 
already been explored, for instance, in [23, 31]. The role of continued fractions in 
this setting is to pick out the correct solutions of the (generalized) Sturm-Liouville 
equations [63], whose relationship to quasidiffusions in general is well-laid out in 
[51]. 

Suppose that Z is a bilateral birth-and-death process. That is, Z is a continuous 
time Markov chain on the integers Z that only makes ±1 jumps. We assume for 
concreteness that Z is killed if it reaches ±oo in finite time, although this assumption 
does not feature in the recurrences we derive in this section. Write (3 n (resp. 5 n ) for 
the rate of jumping to state n + 1 (resp. n — 1) from state n. 

For n G Z, let r n = inf{£ > : Zt = n} be the hitting time of n, with the usual 
convention that the infimum of the empty set is +00. Set 





:= E n [e- XT 




•=E"[e- Ar ' 




:= E"[e- Ar 



Note that 

H n , m (X) = Hl(X)Hl +1 (X) ■ ■ ■ ffl_ 1 (A), m > n, 

and 

H n , m (\) = H^H^X) ■ ■ ■ Hi +l (X), m < n, 

and so the fundamental objects to consider are H* and H^. 

Conditioning on the direction of the first jump, we get the recurrence 

Hi(X) =E" [e-^l Tn _ 1<Tn+1 +e- x ^l Tn+1<Tn _ 1 E n+1 [e^- 1 ]] , 

= 1 — 1 2 1 \ + x — 1 a -l \ H i+iW H n{X), 

On + Pn + X d n + fi„ + A 
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which, putting p n := 4 s -, can be rearranged as a pair of recurrences 
(3.1) H,[(X) 



(3.2) , l (X) = l + Pn + — 



A p n 



This leads to two families of terminating continued fractions that connect the 
Laplace transforms for different values of n, namely, 

Kw - On. 



H n {\) 



1 + Pn 



and 

Hi(X) = l+ P n-l 



Pn-l 



By exchanging S n and (3 n we get similar relations for H^, 



1 + Pn + ±~PnK-lW 
PnHl_ 1 (X) = l+p n + ^-- 



0n Hl{X) 



If we define 



s n (z) := and s n (z) - 



l + p n + ^+Z l+/5„ 1 + ^- + 2 ' 

then we can write the resulting four continued fraction recurrences as 
(3-3) -#i(A) = s„ o s n+1 o ■ • • o s n+m -i(-H* +m (X)), 

(3.4) = sVl o S n _ 2 0...0 ^(--J—), 

1 1 , 

(3.5) - , = S„+i o s„ +2 o • • • o S n+m (-— t — ), 

tin(X) H n+m (X) 

(3.6) -Hl(X) = S„ o s„_x o ■ • ■ o s n -m+i(—Hl_ m (X)). 

In the context of a unilateral birth-and-death chain (that is, the analogue of 
our process Z on the state space N), the context considered in [23, 31], there is 
theory giving conditions under which such continued fractions converge and their 
classical values give the corresponding Laplace transform. In the bilateral case, not 
all of the above continued fraction expansions can converge to the classical values, 
because that would imply, for instance, that H^X) = (Hj l _ 1 (A)) -1 , but two Laplace 
transforms of sub-probability measures can only be the reciprocals of each other if 
both are identically 1, which is certainly not the case here. 

In the next section we consider bilateral chains arising from instances of our 
process £ on T and discuss circumstances in which Laplace transforms of hitting 
times arc indeed given by their putative continued fraction representations. 
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4. Hitting times on a scattered subset of T 

Suppose in this section that for some a G T the infinite set T n (a, +00) is discrete 
with a as an accumulation point. Write Tn (a, +00) = {t n : n G Z} with t n < t n+ \ 
for all n£Z, and define Z : Tn (a, +00) — > Z by Z(t n ) := n. Then the image under 
Z of £ killed when it exits (a, +00) is a bilateral birth-and-death process that can 
"reach —00 in finite time and be killed there". 
From Proposition 2.2, the jump rates of Z are 

1 1 
(4.1) S n = — and /?„ = 



(tn ~ tn-l)(~tn+l ~ tn— l) (t n + l ~ t n )(t n +l ~ t n —l) 

(t n +l ~ tn) 



(4.2) and so p n = 



(tn tn — l) 



The convergence properties of the continued fraction expansions given in (3.3)- 
(3.6) can sometimes be determined by the behavior of Tn (a, +00) in the neighbor- 
hood of its endpoint a. We refer the reader to Section 12 for a review of the theory 
of limit-periodic continued fractions that we use. 

It is clear from the construction of £ as a time change of Brownian motion that 
ia£{t > : & = a} = 0, P a -a.s. and inf{< > : & ^ a} = 0, P a -a.s. That is, a is a 
regular instantaneous point for £. Thus 

(4.3) lim p- fl [ e - AT -"- 1 ] = 1 

n — >oo 

and 

(4.4) lim p-nf e -AT_„ + i J T _ n+i < oo] = 1. 

n — >oo 

Note that /3_ n — > 00 as n — > 00. Suppose further that p_ n — > p G (l,oo) as 
n — > 00, Then s„_ m — ► s* as to — > 00, where s*(z) := 1 _ | ^ ) < i 1+z , a transformation 
with attractive fixed point — p^ 1 and repulsive fixed point —1. It follows from (4.3) 
that lim n _, 00 Jj£ n (A) = 1. So, by Theorem 12.1, -ff£ n is not equal to the classical 
value of the non-terminating continued fraction corresponding to (3.4). Also, by 
(4.3), 

fft n (A) : =E-"[e- AT -+i] 

= P- n {r_„ +1 < oo}E-"[ e - Ar -+ 1 I r_„ +1 < 00] 

= ^-""^^ -"[e-^-^lr-n+i < 00] 
\t-n+l — a) 

1 

— > — as n — » 00. 
P 

Thus, Theorem 12.1, applied with indices reversed, implies that the continued 
fraction expansion in (3.6) converges to the classical value. So, for each n G Z, 
ii^ = — lim m ^oo s„ o • • • o s n _ m (0) = Un/Un+i, where {Uk} is the minimal solution 
in the negative direction to 

(4.5) U k - 1 = (l + pf + hu k -p^ 1 U k+1 . 

Ok 
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In particular, this says that the Laplace transform of the upwards hitting times for 
the process killed at a are given by a simple formula in terms of the {U m }, 

m-l f T 
H n , n+m (\)=]]Hl +k (X) = ^±±, m>0. 



k=0 



Un+1 



Suppose now that (3 n and S n converge to as n — > oo in such a way that p n — > p G 
(l,oo). An equivalence transformation of continued fractions relate the continued 
fraction implied by the recurrence (3.3) and the continued fraction implied by the 
equivalent recurrence 

(4.6) -Pn-iHk{\) 



l3 n + 5 n + \- l3 n H n+1 {X) 

Since f3 n and 5 n tend to zero as n — * oo, the limiting transformation is singular and 
the fixed points tend to zero and —A. Since < H 7 [(X) < 1, lim^oo (3 n -iH n (X) = 
for all A > 0, which is the attractive fixed point of the transformation. Theorem 
12.1 implies that the continued fraction converges to the classical value, which is 
given by the ratio of the minimal solution in the positive direction of the recurrence 

(4.7) V k+1 = (J3 k + 4 + X)V k - 0k-i5 k Vk-i. 

However, Uq := Vq and 



defines a one-one correspondence between solutions to (4.7) and solutions to 

(4.8) U k+1 = (l + p k + ^-)U k - PkUk-L 

Pk 

Note the sequence U is not the same as in (4.5). Since this correspondence maps 
minimal solutions to minimal solutions, if we denote by {Uk} the minimal solution 
in the positive direction to (4.8), then 

tf,i(A) = # l 



U n -i 
and 



m-l ~ 

H„ t n- m (X) = [] H l n _ k (X) = m > 0. 



5. Introducing the process on T q 

Note: For the remainder of the paper, we restrict attention to the state 
space T = T, := {q n : «eZ}U {-q n : n G Z} U {0} for some q > 1. 

In this case the measure p defining the time change that produces £ from Brow- 
nian motion is given by p — p q , where p q ({q"}) = (q n+1 — q"~ 1 )/2, p q {{~ q n }) = 
p q ({q n }), and p q ({0}) = 0. Let | denote the Markov process on T q n (0, oo) = {q k : 
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k £ Z} with distribution starting at x which is that of £ started at x and killed 
when it first reaches 0. 

By Proposition 2.2 the generator Q of £ is defined for all / £ Co(T q ) for which 
the following is well-defined and defines a function in Co(T g ), 



(5.1) (gf)(x) := 



where c q := q^ 1 (q — 1) 2 (1 + <?). In particular, when our process is at any point 
x 7^ 0, it waits for an exponential time with rate proportional to x~ 2 and then 
jumps further from with probability 1/(1 + q) or closer to with probability 

«/(! + ?)■ 

We first reinforce our claim that the process £ on T q is a reasonable g-analogue of 
Brownian motion by showing that £ converges to Brownian motion as the parameter 
q goes to 1. 

Proposition 5.1. For each q let x q £ T q be such that q J. 1. Then the 

distribution of £ started at x q converges as q J. 1 ( with respect to the usual Skorohod 
topology on the space of real-valued cadldg paths) to the distribution of Brownian 
motion started at x. 

Proof. Let (-B*)teR + be a standard Brownian motion with £?o — and let denote 
the jointly continuous local time process of B. 
Set 



JR 

and 

Of := hfi> : A? > t}. 

Then the process (x q + B{9^ ))tgR + has the distribution of £ under F Xq . 

Since /x 9 converges vaguely to the Lebesgue measure m on M as q I 1, we have 



lim< = / £ a - x m(da) =u 

uniformly on compact intervals almost surely, and hence 

limflf = t 

uniformly on compact intervals almost surely. Thus x q + -B((9^ ) converges to x + 
uniformly on compact intervals (and hence in the Skorohod topology) almost 
surely. □ 

The following lemma shows that £ obeys a scaling property similar to that of 
Brownian motion. 

Lemma 5.1. The distribution of the process (£t)tgR + under P x is the same as that 
°f (^£(j 2 t)teR+ under F qx . A similar result holds for the killed process £. 

Proof. The claim for the process £ is immediate by checking that properties (P), 
(II) and (III) hold for (^ q 2 t )teR+- Alternatively, one can verify that the generators 
of the two processes agree, or use the time-change construction of £ from Brownian 
motion and the scaling properties of Brownian motion. The claim for the killed 
process follows immediately. □ 
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Perhaps the easiest things to calculate about the distribution of £ are the mo- 
ments of £t. Formally applying the formula for the generator of £ from Proposition 
2.2 to the function f(x) = x k gives Qf(x) = ^-(1 - q k ){l - q k - 1 )x k -' 2 . As for 
the particular cases of k = 1,2 considered in the proof of Proposition 2.1, we can 
use Dynkin's formula (2.6) and an approximation argument to get the recursion 
formula 

W[$]=x k + f W{Gx k ){^ s )ds = x k + f ^— — (1 — q k ){l - q k ~ l )E x [£ k ~ 2 ] ds, 
Jo Jo c q 

and hence, using the notation introduced in Section 11, 



2|(fc-m) 

where we mean that the sum goes over all < 7« < k with the same parity as k. 
The formula shows that, say for x = 0, the fc th (even) moments grow like 

q i H 1 +°( 1 ))t! . This rate of growth is too fast to guarantee that the moments charac- 
terize the distribution of £t. Note that some well known distributions have moments 
with this rate of growth, for example, the standard log-normal distribution has k th 

k 2 -k 1 

moment e~, as does the discrete measure which assigns mass proportional to e~ 5_ 
at the points e k , k € Z, [12], 2.3e. 

From Proposition 5.1 we would expect informally that the moments of £t should 
converge to those of a Brownian motion at time t as q j 1. Recall that c q = 
q~ l {q — 1) 2 (1 + q) and observe that \\m. q i\{q - l)~ l (q\q)i = (— l) e £l. Therefore, if 
we take x q € T q with lim 9 |i x q = x £ R, then we have 

^ ri i., fc! /1\ 2 (fc — m)! fc-m „, 
limEM^ = > — — - Vi M~x m . 

m— v / \ / V 2 / 

2|(fc-m) 

We recognize the expression on the right hand side as being indeed the k th moment 
of a Gaussian random variable with mean x and variance t. 



6. Hitting time distributions for T q 

We once again stress that for the remainder of the paper we are considering the 
process £ on the state space T q . 

The general considerations of Section 4 apply to T q n (0, oo). In the notation of 
that section, t n = q n for n £ Z. The death and birth rates for the corresponding 

— 2n + l — 2n 

bilateral birth- and-death process on Z are, respectively, — - — and — c — , where we 

recall that c q = q~ x {q — 1) 2 (1 + q). 

To avoid the constant appearance of factors of c q in our results, rather than work 
with £ and its counterpart £ killed at 0, we will work with the linearly time-changed 
processes X = £(c 9 -) and X = £(c 9 -). Of course, conclusions for X and X can be 
easily translated into conclusions for £ and £. 

The corresponding bilateral birth-and-death process on Z has death and birth 
rates S n = q~ 2n+1 and j3 n = q~ 2n . In the notation of Section 4, p = p n = q and 

/ s zl 

Sn[Z) ~ (l + q)+\q^+z- 
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Moreover, we have t„ = inf{i € R+ : X t = q n }. Note, by the scaling properties in 
Lemma 5.1, that Hq(A) = i^(g" 2n A), and # T (A) = Hl(q- 2n X). Moreover, recall 

# n ,„_ m (A) = iJ^(A)JI„_ 1 (A) ■ (A) 



that 



and 

#„, ra+m (A) = ff,,i(A)^ T l+1 (A) • ■ • i/.I +m _ 1 (A), 

so to compute _ff„.„_ m (A) and i?„ jn+m it suffices to compute _ffp(A) and 7Jq(A). 
From Section 4 we have 

(6.1) H^X) = Q - , 

q 

1 + 9 + A 



l + q + Xq 2 - 



(6.2) if J (A) 



1 + q + A</ 4 
1 



o ' 

l + <7 + A 



1 + q + Xq- 



l + q + Xq~ 4 



Closed-form expressions for continued fractions of this form are listed in Ra- 
manujan's "lost" notebook (see the discussion in [5]), and evaluations for various 
ranges of the parameters (although not all the values we need) can be found in 
[5, 32, 33] (although in the last several parameter restrictions are omitted). 

Theorem 6.1. (i) The Laplace transform of the time to go from 1 to q -1 for both 
X and X is 

o( } ~ a oM-Artjj*)' 

An alternative expression is 

01 ' (A^ + i) iM0;-^r; <r 2 ;-^)" 



(ii) The Laplace transform of the time to go from 1 to q for X is 
#o T (A) = 



l KM0;-A<r 3 ;<r 2 ;<T 3 ) 



(g + A) i^i^-Ag- 1 ;?- 2 ^- 3 )' 

Proof, (i) Consider the first expression. Since the continued fraction (6.1) converges, 
by Lemma 12.1 and equation (4.6), the value of q~ 2n -ff^(A) is given by the ratio of 
consecutive terms of the minimal solution to 

W n+1 = ((1 + q)q- 2n + X)W n - q- in+3 W n - V 

This recurrence is found in [32] (but with their q as our q^ 1 ), and the minimal 
solution is shown to be 

U n (X) := q -Mn-D Q_y oM _. . q -i._J_y 
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For the second expression, we evaluate (6.1) as follows. Set 
(6.3) r„(A):= kMO; -^T 2n-1 ; <T 2 ; -jq- 2n ~ 2 )e q -, (-Aq 2 "" 1 ), 



ft re (A) := -<? 



-2n+2 r n(A) 



r„_i(A) 



(6 ' 4) g ^+3 101 ( O; _l g -2n-l ;g -2._l g - 2 n- 2) 

(A + g-3n+l) ^(0; -j-q-2n+l. q -2._l. q -2 n) 

Then, from equation (17) in [5] , 

-4n+3 

(6.5) h, 



(1 + q)q~ 2n + A + h n+ i ' 

This transformation tends to a singular transformation as n — > oo and the fixed 
points tend to x = and y = —A. By (6.4), q 4n h n — > as n — > oo, so /i„ — * 0, and 
convergence to the classical value holds. However, the continued fraction coming 
from (6.5) is related by an equivalence transformation to the continued fraction 
coming from the relation 

q 2(n ^h n = 



l + q + \q 2n + q 2n h n+1 ' 
This is exactly what is needed to evaluate (6.1), and hence 

fli(A) = -q 2 ^K{\) = 

r„-i(A) 

Note that this also shows that 9 - "( n-1 V ra (A) is a minimal solution in the positive 
direction to the recurrence 

U n+1 (X) = ((1 + q)q- 2n + \)U n (X) - g- 4 " +3 C/„_i(A), 

and is hence equal to U n above up to a constant multiple, 
(ii) Define 

r' n {\) = q - n l <t> 1 (0;-M 2n - 3 ;q- 2 ;q- 3 )/e q -*(-M 2n -% 



9nW 



r 



-ra+1 



\i + Aq- 2 "- 1 )- lM °> - A 9" 2n " 1 ;9~ a ;«" 3 ) 



9 H ' iM0;-M- 2n - 3) ;q- 2 ;q- 3 )' 

Equation (13) in [5] simplifies to 
(6.6) g n (X) 



l + q + \q- 2n +g n+1 {\)' 
The fixed points of the limiting transformation arc —1 and —q, and 

1 kMOjO;?" 2 ;^ 1 ) 1 



lim g n (X) 



9 if (o^;?- 2 ;?" 1 ) q 



Hence, by Theorem 12.1, <? n (A) is equal to the classical value of the continued 
fraction implied by (6.6), and r' n is a minimal solution in the negative direction to 
the recursion 

U n+1 = (l + -+\q 2n - 1 )U n --U n - 1 . □ 

q q 
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Let T-oo denote the death time of X. Equivalently, is the first hitting time 
of by X. Write := E 9 " [e- AT —j. 

Corollary 6.1. TTie Laplace transforms of various hitting times for X are given 
by 

q m 2 -2mn „ ^ (_ ; 0; a" 1 ; ) 



1 101 (Pi - Aq2"+i ' » ~ Aq 2 "+ 2 ) 

(-Ag 2 »-l ;(Z - 2 ) m 1 ^ 1 (0-- T -^ Tq 2 m . q -2.__^ q 2 m) 

q 2,— )e q -<-M 2n 



F„,-oo(A) = 1 1 (O;- T ^; (Z -;--^)e 9 - 2 (-Aq-- i )/e ? - 2 (-), 



and 

1 x^^-Ao 2 "- 3 ;?- 2 ;?- 3 ) 



(A) 



Proof. The only result that requires proof is that for H n , t - 00 (X). However, by (11.2) 

lim 1 1 (O;-I O " 2 "- 1 ;q- 2 ;-iq- 2 "- 2 )= ^(0; -; q' 2 ; -) = e q -, {-). 
n— >oo A A q Q 

□ 

We can apply known identities to obtain alternatives to the expressions for the 
Laplace transforms in Theorem 6.1 and Corollary 6.1. For example, equation (13) 
in [51 gives 



r m _ _ i0i(Q;-g *A;g 2 ;g *) 
" i0i(O;-«-iA;(?-2;«-8)- 



#d(A) = l 



Similarly, equation (17) in [5] gives 



The relation 

{w;q)oc i<f>i(0;w;q;c) = (qq)^ i</>i(0; c; g;to) 

follows from (III.l) in [30] upon sending 6^0, letting a = w/z, and sending z — > 0. 
Similarly, the recurrence 



<M0;-A<? fe - 4 ;<r 2 ;<r 3 ) 



r i 



Xq k-^;q j oc fc _ 3 2 . -i 



-2-1 



0i(O;-A fl fc - 3 ;o- 



\Q X ;9 2 )oo 1 _ 2 1 

>i(0;-T-fc53 ;-T-irr) 



comes from (III. 31) in [30] by sending b — ► 0, letting a = w/z, and sending a — > 0. 
Both of these identities can be used to obtain alternative formulae for Hq and i/J . 
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We can invert the Laplace transform flo,— oo m Corollary 6.1 to obtain the dis- 
tribution of the time r_oo for X to hit starting from q n . Note first of all that 



1 



(-Ag 2 "-i;<r 2 ) c 



J! 1 + W^q- 21 



i=0 



n 

i=0 



q 2i-2n+l + X' 



Similarly, 



i^i(0 



1 -.-2. 1 



n+1 ' q ' Ag 2 ™+ 2 ^ 

E [(A<z 2 " +2 ) V^ 15 (-j^T ; ?- a )*(ff- a ; <T 2 )* 



n -1 



fc=0 



fc-1 



En 

fc=o ;=o 



-2(/+n)-l 



(g-2(i+r 



A) 



,-fe 



(<r 2 ;<rV 



Thus under P q the killing time t_oo has the same distribution as the random 
variable 



N 



<f 



iV 

g 2 - E i 2t ~ lT * 



i=l 



where the Ti are independent rate 1 exponentials and N is distributed according 
to a g-analogue of the Poisson distribution [42] , namely, 



P{7V = fc} = 



1 



,-fe 



k > 0. 



It follows that under P 9 the distribution of T-oo is also that of the random 
variable 



i=0 



A partial fraction expansion of the Laplace transform shows that a convolution of 
exponential distributions, where the i th has rate o^, has density 



Hence Y^T=o 1 2j '^j nas density 



•nV ^ 1 ft ( - 1 

11 \l - q -2(k-j) J 11 y 1 _ q - 2 (k-j)J 



(6.7) 



(<r 2 ;g~ 2 )oc ^ (<? 2 ;<? 2 )j 



E 



J'=0 



v(-«- a )E 



(_l)j g -j(j-i) e -<? 2 

(<r 2 ;<r 2 ); 



We note in passing that the random variable Sjlo 9 2j ^i nas ^ ne same distribution 
as the exponential functional of the Poisson process I q -i investigated in [4] (see also 
[3])- 

The following result is now immediate. 
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Proposition 6.1. Under ¥ q " , the hitting time of for X has density 
~JTT £ , f"\, g 2 ("+">+V(*g 2(n+m)+1 ), t > o, 

where f(t) is defined in (6.7). 

Recall that for Brownian motion started at 1, the hitting time of has the 
stable(i) density 

1 ( 1 \ 

: exp — — , t > 0. 



It follows from Proposition 5.1 that the distribution of c q q 2N YHtLa 1~ 2l Ti converges 
to this stable distribution as q J, 1. From Lai's strong law of large numbers for 
Abelian summation [53] we have that 



limE(! - q~ 2 )q~ 2l T t = E[T ] = 1, a.s. 



<tl- 

and so c q (l — q~ 2 )~ 1 q 2N also converges to the same stable distribution. Taking 
logarithms, we obtain the following result. 

Proposition 6.2. As q J, 1, the distribution of the random variable 

2(lo gg )JV + Iog(g-l) 
converges to the distribution with density 

cxp ( — -(x + exp(— x)) ) , — oo < x < oo. 



7. Excursion theory for Tg 

Recall that X — ^(c q -) under P x has the same distribution as (B(0 C t))teM+, where 
B is a Brownian motion started at x with local time process £, 9 is the right- 
continuous inverse of the continuous additive functional A t = J £f (j,(da), and fi q 
is the measure supported on T q that is defined by H q ({q n }) = (q n+1 — q n ~ 1 )/2, 
fi"({-q n }) = fi(U n }), and ^({0}) = 0. 

Recall also that is a regular instantaneous point for X. Thus X has a continuous 
local time L at that is unique up to constant multiples. We can (and will) take 
L t = £g t . The inverse of the local time is a subordinator (that is, an increasing 

Levy process). Also, there is a corresponding Ito decomposition with respect to the 
local time of the path of X into a Poisson process of excursions from 0. In this 
section we determine both the distribution of the subordinator (by giving its Levy 
exponent) and the intensity measure of the Poisson process of excursions. 
We begin with the following result, which is immediate from Lemma 2.1. 

Lemma 7.1. The process X is reversible with respect to the measure [i q . In par- 
ticular, fi q is a stationary measure for X . 

We use the excursion theory set-up described in Section VI. 8 of [61], which we 
now briefly review to fix notation. Adjoin an extra cemetery state d to T q . An 
excursion from is a cadlag function / : R + — > T q U {d} such that /(0) = and 
/(*) = d for t > C, where C := inf{* > : f(t) = d or f(t-) = 0} > 0. Write U for 
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the space of excursion paths from 0. Using the local time L, we can decompose that 
paths of X under P x into a Poisson point process on M. + x U with intensity measure 
of the form m <g> n, where m is Lebesgue measure and n is a cr-finite measure on 
U called the ltd excursion measure. The measure n is time-homogeneous Markov 
with transition dynamics those of X killed on hitting (and then being sent to d). 
Thus n is completely described by the family of entrance laws n t , t > 0, where 

n t (r) := n({/ G U : f(t) G T}), T C T,. 

Let R\ denote the X-resolvent of X for A > 0. That is, 

~ M V x {X t G T}dt 



Jo 



for x G T q and rcT,. In order to identify n. we begin with the following general 
excursion theory identity (see equation (50.3) in Section VI. 8 of [61]). 



(7.1) 

where 

(7.2) 



-XI 



n t ({y})dt = R x (0,{y}) 



-As 



dL s 



Now, setting T x := inf{f G R+ : X t = x}, x G T 9 , 

}) 

A^(M) 



i? A (0,M)= hmi? A (x,M) 

.r — 



= iim n r .)ri^ A ^w) 

■t^O ^({.t}) 



lim 



/i 9 ({x}) 

hm 

^({2/})^ \e- XTo ] lim flA(0 ; { ^ }) 



E y 



R\{x, {x}) 



//'({j/})— E» [e- ATo l liniE 



^({j/})!" [e~ AT °l -E° 
M 9 (M)E y [ e - ATo ] E° 



e~~ As d£* 



- xt de° 



d£ 



where we used Lemma 7.1 in the second and sixth lines, and a change of variable 
in the final line. 



Thus, 



so that 



-A/ 



n t {{y})dt = pP{{y})W [e~ XT °] 



nt({y})=» q ({y})- 



{To € dt} 
dt 



6.1. 



Now E«" [e 



ATo 



] = fl n) _ 00 (A), and so we obtain the following from Corollary 
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Proposition 7.1. The family of entrance laws (n t ) t> o is characterized by 

poo 

/ e- xt n t ({q n })dt 
Jo 



= ^ 1 (9 2 -l)i0i(O 



1 .-2. 1 



\ q 2n+i'1 ' Ag 2„+2 



and n t ({-q n }) = n t ({q n }), n g Z. 

Let 7 denote the right-continuous inverse of the local time L, so that 7 is a 
subordinator. Thus E°[e~ A7t ] = e~*^ A ) for some Laplace exponent ip. 

Proposition 7.2. The distribution of the subordinator 7 is characterized by 
A (q 2 - 1) (-1, -A.T 2 ; g- 2 )^ _X(q 2 - 1) e,- a (-£)e ff - 3 (-£) 



^(A) 



9(-|,-^;?- 2 )= 



q ' Aq 

Proof. We note the relationship 

k a = E° 



Q e q- 2 (-\) e q- 2 (-M 2 ) 



dt 



dt 



Hence, from equation (7.1), 



^(A) = A/ e- M n t {T q \{0})dt = 2X 



2A^>- 1 (g 2 -l) 1 1 (O;-^ 



e- At n t (T 9 n (0,oo)) (it 
1 



<T 2 ; 



A? 



2n+2 



I - 



Using the following identity to simplify the sum, 

(-Ag 2 "- 1 ;?- 2 )^ = (--^U-Ag 2 "" 1 ;?- 2 ^ 

= (--;g- 2 )oo(-4;^ 2 ) Il 9" 2 A", 

9 ?A 

we can write part of the above as 

g n i^ilO;-!^;?' 2 ;-!^) 



E 

n6Z 



(_Ag2n-l. g -2) 
1 



^ q — n{ri — 1) ^ — n 

— (° 



1 



1 



( g ;s 2 )oo nez ( q \'i<i 2 ) 



1 



A q2n+1 '3 ' ' Ag 2n_t 
fc(fe-l)-fe(2n+2) ^-n-fe 



n£Z,fc>0 gA' y / n V Ag 



H;<r 2 ) 



i-^i- 2 )n{- q -rr;q- 2 Uq- 2 ;q- 2 )k 
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and by changing indices we get 



q— n(n— 1)— fc(fc-l)— fc(2n+2)^\ — n-k 
Jt>0 ("i; 1- 2 U-^; q- 2 )k(q- 2 ; <T 2 ) fe 

~ — mini— 1) \ — m rl —2k 

E g 'A \ ~» 1 

(_ 1 • n -2\ 2-^1 ' 



= o^i(-;--^]q~ 2 -,-j)e q -2{q~ 2 ). 

Recall that e q -i{z) = l/(z;q^ 2 ) oc . Moreover, using equation (11.3), we can 
rewrite the oV'i as a product. 

l.__ 2 . 1, (q- 2 ,-{,-\q- 2 ; q- 2 )oo 



oipi(-;-T-'9 ;-t) 



V y ' ^ (-^.J;«- a )co 

The result now follows. □ 

It follows from the scaling property Lemma 5.1 and the uniqueness of the local 
time at up to a constant multiple that (L 9 2 t ) tg R + has the same distribution 
under P° as a constant multiple of L. Consequently, the exponent ip must satisfy 
the scaling relation i/j(q~ 2 X) = cip{\) for some constant c. Note from the formula 
in Proposition 7.2 that, indeed, 



8. Resolvent of the killed process on T q n (0, oo) 



Let R\ denote the resolvent of the process X on T q n (0, oo) killed at 0. Recall 
that X goes from q n to at rate q- 2 ^ 1 and from q n to q n+1 at rate q~ 2n . 

Thus the exit time from q n is exponentially distributed with rate q~ 2n+1 + q~ 2r \ 
the probability of exiting to q n ~ x is and the probability of exiting to q n+1 is 

Moreover, E«" _1 [ e - AT "] = ^_ X (A) and E«" +1 [ e - AT "] = i^ +1 (A). From the 
strong Markov property we get the recurrence 

Rx(q n ,{q n }) 

f Rx(q n ,{q n }) 



A 



1 

-2n+l 



5 



A + fy -2»+l 

so that 

Rx{q n ,{q n }) 

= {\+(q- 2n+1 



q 2n q 



1 



-2n+l 



\ + q 



-2n+l 



q 2n q 



-Hi +1 (X)), 



') 



1 - 



Substitute any of the explicit formulae for Hj l _ 1 and H^ +1 from Section 6 to get 
an expression for the on-diagonal terms of the resolvent in terms of basic hypcrgc- 
omctric functions. 
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To obtain the off-diagonal terms, use the observation 

R x (q m , {q n }) = E« m [e" AI >] R x (q n , {q n }) = H m>n {X)R X {q n , {q n }) 

and then substitute in explicit formulae for -ff m>ra (A) from Section 6 to get expres- 
sions in terms of basic hypergeometric functions. 

Ideally, one would like to invert the Laplace transform implicit in the resolvent 
to obtain expressions for the transition probabilities P x {Xt = y}- We have not been 
able to do this. 

9. Resolvent for T q 

Recall that R x is the resolvent of the process X. The resolvents R\ and R x are 
related by the equations 

f R x (x,{y})+E x [e" AT °] R X (0,{y}), x,y € R+, 
Rx(x, {y}) = I W [e- XT °] R x (0, {y}), x>0,y<0, 
{R x (-x,{-y}), x<0. 

Recall equation (7.1), which says that R x (0,{y}) = k x f °° e~ xt n t ({y}) dt. We 
know from the proof of Proposition 7.2 that k x = and the statement of Propo- 
sition 7.2 gives a simple expression for ?/^(A) as a ratio of infinite products. Proposi- 
tion 7.1 gives an expression for / °° e~ xt n t ({y}) dt in terms of basic hypergeometric 
functions. Again noting that E 9 ™ [e~ AT "' 1 ] = H m , n (X), we substitute in explicit 
formulae for H m n {\) from Section 6 to get expressions for R x (x,{y}) in terms of 
basic hypergeometric functions. 

10. A remark on spectral representations 

An alternative approach to finding explicit formulae for the quantities of interest 
would be to find a spectral representation for the generator. This is well-described 
for general quasidiffusions by Kiichler and Salminen [51], who build on the spectral 
theory of strings [13, 41, 63]. Once one has found solutions to the Sturm-Liouville 
equation Gu = — Xu with appropriate boundary conditions, and the orthogonalizing 
(spectral) measure, one can write down explicit formulae. 

One possible method for carrying this out is to use the well-known spectral 
representation of transition probabilities of a unilateral birth-and-death process, 
for which the appropriate cigcnfunctions arc a family of orthogonal polynomials 
(see, for example, [37, 38, 39, 62]). If we kill X at q~ n for some n G Z to obtain a 
process on {q~ n+1 , q~ n+2 , . . .}, then the corresponding unilateral birth-and-death 
process has a specialization of the associated continuous dual q-Hahn polynomials 
as its related family of orthogonal polynomials [32] . However, we have not been able 
to "take limits as n — > oo" in the resulting spectral representation of the transition 
probabilities to obtain similar formulae for X. 

Note that our expression for the density of the hitting time to zero of Proposition 
6.1 appears to be close to a spectral decomposition - compare to Theorem 3.1 in 
[51], which gives the density as 



- [ e- x2t C(x;X)p(d\) 
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where C is a particular solution to the Sturm-Liouville equation and p is the spectral 
measure. However, to put our expression in this form, the two summations need to 
be exchanged, which is not straightforward. 

11. Background on basic hypergeometric functions 

For the sake of completeness and to establish notation, we review some of the facts 
we need about basic hypergeometric functions (otherwise known as q -hypergeometric 
functions). For a good tutorial, see the article [45] or the books [2, 30]. In order to 
make the notation in our review coincide with what is common in the literature, 
take < q < 1 in this section (this q usually corresponds to q~ 2 in the rest of 
the paper). 

Define the q-shifted factorial by 

n-l 

(z; q) n := Y[ (1 - zq k ) for n G N, Z £ C, 

k=a 

oo 

(z; q)oo ■= n ( x Z( i k ) for \ z \ < L 

k=0 



The definition of (z; q) n may be extended consistently by setting 
(z;g) fc = , {z \ q)o ° forfceZ, zeC. 

(zq k -,q) oa 

It will be convenient to use the notation 



(oi,o 2 , • ■ .,a r ;q) k = {a 1 ;q) k (a 2 \q)k ■ ■ ■ (a r ;q)k- 

The q-hypergeometric series are indexed by nonnegative integers r and s, and 
for any {a^ C C, {bj} C C \ {q~ k }k>o are defined by the scries 

, , , ^ (a 1 ,... 1 a r ;q) k ((-l) k g h ^y+^ z k 

r<p s {ai,...,a r ;bi,...,b s ;q;z) := > j- . 

^ (bi,...,b s ,q;q) k 

Note the factor (q;q)k on the bottom, which is not present in the definition used 
by some authors. The series converges for all z if r < s, on |z| < 1 if r = s + 1, and 
only at z = if r > s + 1. Using the property that 

lim^ = (-l)V^, 

a — !-oo d 

we get the following useful limit relationships 

z 

(11.1) lim r+ i0 s (a, oi, . . . , a r ; h,..., b s ; q; -) = r <j> a (ai, . . . , a r ; b u . . . , b 3 ; q; z), 

a^oo a 

(11.2) lim r (p s+1 (a 1 , . . . , a r ; b, b 1} . . . , b s ; q; bz) = r (j) s (a 1 , . . . , a r ; h, . . . , b s ;q; z), 

b— *oo 

as long as the limits stay within the range on which the series converge. 
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(az;q) c 
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i<Ma;-;g; z) 



if\z\ < l,\q\ < l,oG 



There are (at least) two commonly used g-analogues of the exponential function. 



e q {z) := !0 o (O;-;9;z) = 1 = ^ 



and 



:= o^o (-; -;<?; -2) = 



e <?(- z ) 



(?;?)« 



(-z;g) c 



for I2I < 1, 



fc=0 



fe(fc-l)/2 



(-*)* 



(?;<?)* 



for z G 



The bilateral q-hypergeometric series also appear in our results. They are defined 



by 



ri) 3 (ai > ... ) a r ;b 1) ...,b a ;q;z):= 

fc= — OO 



(g 1 ,...,a r ;q) k ((-l) k q^—) s r z k 
(61, . . .,b s ;q) k 



The sum converges for 





•i>8 















< \z\, if s > r, 

< \z\ < 1, if s = r, 



and diverges otherwise. 

We use the following extension of the Jacobi triple product identity (see equation 
(1.49) of [45]) 



(11.3) oi>i{-; c;q;z) := ^ 

k— — 00 



(-l) k q k ( k 1 )/ 2 z k _ (q,z,q/z;q) c 



z > c 



(c;q)k (c,c/z;q)oo 

12. Background on recurrence relations and continued fractions 

For nonzero complex numbers a n and b n , n G Z, consider the three-term recurrence 
relation 

(12.1) U n +i = b n U n - a n U n -i- 

Its connection to continued fractions can be seen immediately by rearranging to get 



U n -1 b r , 



U n+1 



In other words, the sequence W n = U n /U n -i solves the recurrence 
(12.2) W n W n+1 = b n W n - a n . 
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Iterating this recurrence, we get that for any k > 0, 

— a„ 

-W n = 



On 



a n +2 

On+1 — 



u n+2 - ■ . a n+k 

b n+k -W n+k 

We refer to this expression as the continued fraction expansion associated with the 
recurrence (12.2). 

A solution (C/ n ) ng z to (12.1) is said to be a minimal solution if, for all linearly 
independent solutions V n , limn—^ |?7 n |/|V^| = 0. The minimal solution to (12.1), 
if it exists, is unique up to a constant multiple [56]. 

For clarity, define the linear fractional transformations 



Sn{z) = — , 

On + Z 

and write their compositions as = s m +i ° s m +2 ° • • • ° s n and S n = Sq. The 
iassical approximants to the nonterminating continued fraction (sometimes written 
Kf^p 2 -]) are given by 



c 



If we let P„ and Q n be two solutions to (12.1) with initial conditions P_i = 1, 
Pq = 0, Q-i = 0, and Qq = 1, then it is easy to see that 

S n (z) Pn + zPn ~ 1 



Qn ~t~ zQn—1 



The continued fraction Kf^p 2 -] is said to converge if S n (0) converges to a (finite) 
limit as n tends to infinity. If this limit exists, it is called the classical value of 
the continued fraction. However, this is a bit arbitrary, because it can happen, for 
instance, that for all sequences (w n ) n gN that stay away from zero, S n (w n ) converges 
to the same limit, different from the limit of ^(O). The problem is easy to see: 
suppose that a n — > a* and b n — > b* as n — > oo, so that s n — > s*. Each s n has a pair 
of fixed points that converge to the fixed points x and y of s* - suppose |x| < \y\, so 
that x is attractive and y is repulsive. One might imagine that as long as w n stays 
away from the repulsive fixed point of s*, then linin^oo S n \{w n ) must converge to 
x as m — > oo, in which case 

lim S„(w„) = lim S m o 5™ (w„), Mm > 0, so 

n — ^oo n — >oo, n>m 



= lim lim 5 m o5"(w„) 

m—*oo \ n— *oo, n>m 

= lim S m (x) . 

m — >oo 

Note that just by setting w n = S^ 1 (z), we can get S n (w n ) converging to any limit 
in C we'd like - but to do this, the w n we choose must converge to the repulsive 
fixed point. The precise sense in which w n must "stay away" from the repulsive 
fixed point is given in the Theorem 12.1 below. 
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The case in which a„ — > a* and b n — > &*, where if a* = then &* ^ 0, is 
called the limit 1-periodic case. Moreover, if the fixed points of s* are distinct and 
have different moduli the continued fraction is of loxodromic type. All the continued 
fractions we deal with fall into this category. The following combines Theorem 4 
in Chapter II and Theorem 28 in Chapter III of [56]. Here d(-,-) is the spherical 
metric on C. 

Theorem 12.1. Let if[^ 2t ] be limit 1-periodic of loxodromic type. 

(i) There exists an f £ C such that for every sequence (w n ) n gs for which 

liminfn^oo d(w n , >S , ~ 1 (oo)) > 0, when f ^ oo, 
liminfn^oo d(w n , S'~ 1 (0)) > 0, when f = oo, 

we have S n (w n ) — > /. In particular, S n (0) — > /. 

(ii) Let V Q eC and V n = S^CVq) forn > 0. Consider V e C and V n = S' 1 ^) 
for n > 0. Write s* = lim^oo s„ and suppose that s* has fixed points x,y with 
\x\ < \y\, so that x is attractive and y is repulsive for s* . 

• IfVo=f, then linijj^oo V n = x. Moreover, if f ^ oo, then Vq = Uo/U-i, 
where U n is a minimal solution to (12.1). 

• Otherwise, linin^oo V n = y. 

For a proof, see [56]. This implies the following lemma, which we also use to 
introduce some more notation. Note that the relation — W n = S~ 1 (—Wq) for n > 
is exactly the relationship implied by (12.2). Note also that this gives explicitly 
the value of the continued fraction, if it converges to a finite value, in terms of the 
minimal solution to the associated recurrence relation, a result known as Pinchcrlc's 
theorem [56]. 

Lemma 12.1. Suppose that (W n ) n £Z solves (12.2) and that the limits 
(3± := lim \ (b n ± y/b^ - 4a n ) , 

n— too Z V / 

exist, are finite, and the branches of the square root are chosen so that |/3_| < 

If linin^oo W n then lmin—,00 W n = —ft-, and for any fixed m £ Z, the 

sequence U n defined by 

(ni= m+ iW k , n>m, 
U n := s 1, n = m, 

l(nr=n + i^) _1 . n<m, 

is a minimal solution to (12.1). 

Proof. Since /3_ is the limit of the attractive fixed points of the corresponding 
transformations, and f3 + is the repulsive fixed point, Theorem 12.1 says that W n is 
equal to the classical value of the continued fraction 

a„ 
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if and only if lim^oo d(— W n +k> 0+) > 0, so a minimal solution U n exists, and 
W n = ^ n . By definition. U n = j^- for all n£Z. This proves the lemma. □ 

Two continued fractions are said to be related by an equivalence transformation 
if their sequences of approximants are the same. For example, let Cfe, k G Z, be 
nonzero complex numbers. Since for all n > 0, 



c a a Q 



b 



CQ&O 



cociai 



c n b n -\-c n w n 



we say that the continued fraction expansions on either side are related by an 
equivalence transformation. 

Note that since we allow the indices in (12.1) and (12.2) to take values in Z, 
by reversing indices we get another recurrence, another continued fraction, another 
minimal solution, etc. When we need to distinguish, we will refer to, say, U n as 
a minimal solution to (12.1) in the positive direction if the above definition holds, 
and a minimal solution to (12.1) in the negative direction if lim n ^oo U- n /V- n = 
for some (and hence any) other linearly independent solution V n . 
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